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ABSTRACT 

We present submillimetre observations of the / = 3 ^ 2 rotational transition of ^^CO, 
^^CO and C^^O across over 600arcmin^ of the Perseus molecular cloud, undertaken with 
HARP, a new array spectrograph on the James Clerk Maxwell Telescope. The data encompass 
four regions of the cloud, containing the largest clusters of dust continuum condensations: 
NGC1333, IC348, L1448 and L1455. A new procedure to remove striping artefacts from the 
raw HARP data is intro duced. We compare th e maps to those of the dust continuum emission 
mappe d with SCUBA ([Hatchell et al.||2005| ) and the positions of starless and protostellar 



cores ( Hatchell et al.|2007a| ). No straightforward correlation is found between the masses of 



each region derived from the HARP CO and SCUBA data, underlining the care that must be 
exercised when comparing masses of the same object derived from different tracers. From the 
^^CO/C^^O line ratio the relative abundance of the two species ([^^CO]/[C^^O] ^ 7) and their 
opacities (typically T is 0.02-0.22 and 0.15-1.52 for the C^^O and ^^CO gas respectively) 
are calculated. C^^O is optically thin nearly everywhere, increasing in opacity towards star- 
forming cores but not beyond Tig ~ 0.9. Assuming the ^^CO gas is optically thick we compute 
its excitation temperature, Tgx (around 8-30 K), which has little correlation with estimates of 
the dust temperature. 

Key words: submillimetre - stars: formation - ISM: kinematics and dynamics - ISM: indi- 
vidual: Perseus. 



1 INTRODUCTION 

Recent advances in telescope instrumentation have provided an un- 
precedented view of the star formation process inside molecular 
clouds. Near-infrared imaging from e.g. Spitzer (e.g. Eva ns et al.| 
|2009 ) captures the youngest stellar objects, whilst (sub)millimetre 
continuum imaging with e.g. SCUBA or MAMBO maps the very 
earliest stages of star formation, often before a well-defined central 
object is estab lished (e.g. Motte, Andre & Neri|| 1 998 1 1 Johnstone] 
leTa l. 2000, Hatchell et al.|2005||Enocirera l. 2006, Di Francesco 
|et al.^.2008j . Simulations of star formation have also become in- 
creasingly sophisticated (e.g. |Klessen & Burkert||2000[ [Bate &| 
|Bonnell|[2005l |Bate||2009a] with the latest models including tur- 
bulence, magnetic fields and radiative transfer (Price & Bate 2008 , 
|Bate|2009b| . The true, underlying structure in star-forming clouds, 
essential for comparison with models, is hard to disentangle from 
continuum images alone ( [Ballesteros-Paredes & Mac Low||2002{ 
[Smith, Clark & Bo nnell 2008b ): density enhancements along the 
line of sight may superpose and limited spatial resolution and sensi- 
tivity inevitably blends nearby objects together. Spectral lines from 
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molecular species can provide crucial information about the kine- 
matics of and physical conditions inside molecular clouds. They 
may for instance allow the separation of multiple objects, moving 
at distinct velocities, along a line of sight (e.g. |Kirk, Johnstone &| 
|Tafalla|2007| ). Additionally, line emission is arguably the best dis- 
criminator between different support mechanisms for star-forming 
cores, i.e. whether magnetic fields or turbulence dominate the re- 
quired suppression of star formation in the bulk of clouds. 

This paper presents a new survey of the kinematics of star for- 
mation in the Perseus molecular cloud. We observed four fields, 
totalling approximately 600 arcmin^ in the / = 3 ^ 2 rotational 
transition of ^^CO, ^^CO and C^^O using HARP (Heterodyne Ar- 
ray Receiver Programme), a new array spectrograph operating be- 
tween 325 and 375 GHz on the James Clerk Maxwell Telescope 
(JCMT, see [Smith et al.|[2003l [2008a[ [Buckle et all[2009j . These 
data represent some of the largest and highest-quality maps of the 
higher-/ transitions of CO and its isotopologues with good angular 
resolution (< 20arcsec). This paper is the first in a planned series 
which will examine the gas structure and core kinematics (Paper 
11), molecular outflows (Paper 111) and the physical conditions in the 
cloud cores (Paper IV). This paper is organized as follows: the re- 
mainder of the introduction outlines the survey and its aims before 
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Table 1. The transitions of CO and its isotopologues observed. 



Molecule 


A/ 


a 


^trans ^ 

(GHz) 


trans 

(K) 


^crit 

(cm~^) 


CO 


3^ 


2 


345.796 


33.2 


2.54 X 10^ 


13CO 


3^ 


2 


330.558 


31.7 


2.22 X 10^ 


ci^o 


3^ 


2 


329.331 


31.6 


2.19 X 10^ 



^ Transition quantum numbers. 
^ Rest frequency. 

^ Height of the upper / level above ground. 

^ Estimated critical density, assuming a collision cross-section a = 
10~^^ cm~^ and velocity 1 kms^^ so ^crit = A x 10^^ cm~^. The Einstein A 
values were taken from the LAMDA database jSchoier et al.|2005) . 

Sections |2] and |3] detail the observations and data reduction proce- 
dure including an algorithm to calibrate variable sensitivities across 
the HARP array for molecular cloud data. Section|4]presents maps 
of the integrated intensity of all of the isotopologues while some 
initial results including the cloud opacity and excitation temper- 
ature, assuming local thermodynamic equilibrium (LTE) are pre- 
sented in Section[5] We summarize this work in Section[6l 



1.1 The survey 

The transitions of CO and its common isotopologues in the 
345 GHz atmospheric window trace denser and/or warmer regions 
than their lower / relatives - compare the temperature of 5 K and 
density of /icrit ~ 700 cm~^ required to collisionally excite the 
^^CO / = 1 ^ line at 1 15 GHz to Table[l] These conditions are 
similar to those inside star-forming SCUBA cores, allowing us to 
probe the motions of bulk gas in their vicinity at matched resolu- 
tion to SCUBA 850 ^tm observations. This survey is closely linked 
to the JCM T Gould Belt Legacy Survey (GBS, |Ward-Thompsonl 
|et aLpOO"?] ). Our data will eventually be included in that analysis 
and we share similar science goals, namely: 

(i) Map any high- velocity outflowing gas in the I^CO / = 3 ^ 
2 line to investigate mass-loss and accretion in a large sample of 
sources. The presence of an outflow can differentiate between a 
protostellar or starless continuum core, establishing the source age. 

(ii) Establish the level of non-thermal support inside star- 
forming cores from the C^^O transition, which may also be used 
to evaluate the levels of CO depletion. 

(iii) Investigate the turbulent structure of the molecular gas. 

Perseus (see Fig. [TJ is an intermediate star-forming environ- 
ment between low-mass, quiescent Taurus and high-mass, turbu- 
lent Orion (e.g. |Ladd, Lada & Myers|1993] ). The complex consists 
of a series of dark clouds at approximately 3^30"^, +31° with an 
angular extent of about (1.5 x 5)deg, totalling ~ 1.7 x 10^ M0 
( [Bachiller & Cernicharo|p^86ap . The cloud is associated with 
Per 0B2, the second closest OB association to the Sun, with an 
age < 15Myr (see [Bally et al.||2008| ). The association has blown 
a 20deg diameter shell of atomic hydrogen into the interstellar 
medium with the Perseus molecular cloud embedded in its western 
rim. The region is well- studied both as a whole and in individual 
sub- sections with a wealth of a ncillary data availab le, particularly 
from the COMPLETE projecf] ( [Ridge et al.|2006bt . 



Notably, large-area surveys of the dust continuum emission 
in Perseus have been completed with SCUBA at 850 yim ( |Hatchell| 
let al. 12 005 ' "Kirk , Johnstone & Pi Francesco|20 06 ) and Boloc am at 
1.1 m m (Enoch et aL|2006t . The 3 deg^ SCUBA survey ( Hatchell| 
|et al.||2005| ) identified 92 submillimetre cores above their com- 
pleteness limit of 0.4 M0 (in a 14arcsec beam). 80 per cent of 
these cores were grouped in six clusters: NGC1333, IC348, L1448, 
L1455, Bl and B5. The larger and more sensitive Bolocam sur- 
vey (7.5 deg^ down to 0.18 Mp; [Enoch et al. 2006) identified 122 



compact cores. However, only 5-10 cores were in areas not ob- 
served with SCUBA by ' Hatchell et al.| illustrating that most of 
Perseus is devoid of active star formation. In this survey we target 
the clusters of star-forming cores (see Section [2]). Spitzer has pro- 
vided a census of deeply embedded young stellar objects (YSOs) 
in many nearby molecular clouds (e.g. [Evans et al.|2009| ) and this 
population has been associated with continuum cores identified by 
SCUBA (j0rgens'^eFaL|2OO7a 2008) or Bolocam ([Enoch et al" 
2008 2009). In this survey, we use the catalogue of Hatchell et al. 
(2007a), who identified the SCUBA cores from Hatchell et al. 



( 2005) as starless. Class or Class I protostars on the basis of their 
spectral energy distributions (SEDs), which incorporated Spitzer 
fluxes from the IRAC wavebands [3 



Estimates of the distance to Perseus vary from 220 pc ( Cernis 
[1990] ) to 350 pc ( [Herbig & Jones|1983] ). Larger distances are often 
based on the Perseus OB2 association which has an established dis- 
tance of ~ 320 pc from Hipparcos ( d e Zeeuw et al.|199"9] ). However, 
there is some evidence that Per OB 2 may lie behind the molecular 



clouds of interest which are probably at closer to 250 pc (Cernis 
[1993] ). In fact, it may not be appropriate to use a single distance 
for the entire cloud and extinction studies suggest increasing dis- 
tances fro m 220 to 26 pc from west to east (previous references 
and Ce rnis & StraizyS|2003j . Furthermore many authors have sug- 
gested that Perseus is a superposition of a least two smaller clouds 
- the closer is thought to be an extension of Taurus with the more 
distant a shell-like structure (e.g. Ridg e et al.[2006a| ). The latest and 
possibly most reliable measure of the distance to the H2O maser in 
NGC1333 SVS13 found (235 ± 18) pc using its parallax at radio 
frequencies ( [Hirota et al.|2008| ). In this survey we assume Perseus 
to be a single entity at a distance of 250 pc for consistency with the 
majority of recent studies e.g. the Spitzer Cores to Disks team, c2d 
([Evans et al.|2003[[2009| ) and the Bolocam Perseus survey ( [Enoch[ 
[et al.|2006| ). 



2 OBSERVATIONS 

We selected the four largest clusters of continuum cores - 
NGC1333, IC348, L1448 and L1455 - to map in the ^^CO, ^^CO 
and C^^O/ = 3^2 lines (see Fig.[l] the exact areas are detailed in 
Table [2]). All the regions were observed in the three tracers except 
for NGC1333 where only ^^CO and C^^O data were taken since 
it had already been observed in ^^CO / = 3 ^ 2 with HARP (J. 
Swift, personal communication). 

The data were taken as part of the Guaranteed Time pro- 
gramme for the HARP instrument team, over a period of nine nights 
between 17th December 2007 and 12th January 2008 with one map 
of NGC1333 taken on 28th July 2007. The data comprise -43 hrs 
total observing time, approximately 25 hrs on sky translating to an 



^ the coordinated Molecular Probe Line Extinction Thermal Emission 
survey of star-forming regions, see http://www.cfa.harvard.edu/ 
ICOMPLETEI ^ 



^ IRAC (Infrared Array Camera) on Spitzer has four channels at 3.6, 4.5, 
5.8, and 8 i^m. 
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Figure 1. Overview of the molecular gas in the Perseus molecular cloud. The grey-scale is C^^O / = 1 ^ integrated intensity in Kkm s~^ from the FCRAO 
14-m telescope ( Hatchell et al. 2005). The dotted contour encloses regions where the visual extinction is > 3.5, from the map of the COMPLETE team ( RidgeJ 
|et al . 2006b I. The four fields observed in this survey are enclosed in boxes: NGC1333, IC348, L1448 and L1455. Also marked is the Bl molecular ridge. B5 
the other large cluster of SCUBA cores lies roughly 45 arcmin to the northeast. 



Table 2. The regions observed. 



Region 


Centre ^ 


Dimensions 


Area 




RA 


Dec 


Width 


Height 






(h m s) 


(° ' ") 


(arcsec) 


(arcsec) 


(arcmin^) 


NGC1333 


03:28:56 


+31:17:30 


760 


900 


190 


IC348 


03:44:14 


+31:49:49 


1100 


650 


199 


L1448 


03:25:30 


+30:43:45 


800 


500 


111 


L1455 


03:27:27 


+30:14:30 


880 


520 


127 



^ Position of the map centre in J2000 coordinates. 



observing efficiency of 58 per cent. The weather conditions were 
good to excellent throughout with median receiver and system tem- 
peratures of Trx = 180K and Tgys = 325 K for ^^CO compared to 
Trx = 179 K and Tsys = 415 K for ^^COIC^^O. 

The HARP imaging array comprises 16 SIS detectors, ar- 
ranged on a 4 X 4 grid, separated by 30 arcsec. This under-samples 
the focal plane with respect to the Nyquist criterion {X/2D is 6.0 
and 6.2 arcsec at 345 and 330 GHz respectively), therefore standard 
HARP observing strategies take data points in between the nomi- 
nal detector positions to produce a fully-sampled map (see |BuckTel 
|et a l. 2009). The observations for this survey were taken in the 
raster mode at the default sample spacing of 7.3 arcsec. The ^^CO 
and C^^O observations were taken simultaneously, utilizing the ca- 
pability of the new back-end correlator, ACSIS (Auto-Correlation 
Spectral Imaging System, [Buckle et al.|2009 ), to split its bandpass 
into two separate sub-bands. Both sub-bands provide ~250MHz 
of bandwidth on the two lines with a channel spacing of 61 kHz 
corresponding to approximately 0.05 km s~^ The ^^CO observa- 
tions were also taken using a sub-band mode, with both centred at 
the same frequency: (i) for high- velocity gas with a bandwidth of 



~1 GHz, broken into channels 977 kHz wide (0.8kms~^) and (ii) 
to trace the gas in detail with ~250 MHz bandwidth at a channel 
spacing of 61 kHz (0.05 km s~^). 

2.1 Scan strategy 

During raster or "on-the-fly" mapping, the telescope scans contin- 
uously along a set direction (usually parallel to the longest side 
of the map), dumping the data gathered at discrete intervals - as 
little as 100 ms apart with HARP/ACSIS. Once the telescope has 
scanned the length of the map, it steps perpendicular to the scan 
direction and begins a new row. The HARP array is angled at 
arctan(l/4) ~ 14 deg with respect to the scan direction to produce 
a fully- sampled image on a 7.3 arcsec grid ( [Buckle et al. [2009 ). 

The detectors in HARP have varying responses across the ar- 
ray. Additionally, in any one observation the data from a number of 
detectors can be unusable with poor baselines. These in combina- 
tion can lead to maps with variable noise or entirely blank strips. 
Our observing strategy employed three factors to reduce these ef- 
fects: 

(i) The perpendicular spacing between scan rows was set to be 
half the width of the tilted array. Thus, in any row, half of the de- 
tectors are going over previously observed positions and half are 
scanning new parts of the sky. 

(ii) The most effective way to 'fill-in' entirely missing rows is 
to 'basket- weave' i.e. re-observe the map with the scan direction 
perpendicular to the original. 

(iii) A small offset between map repeats, of one pixel's length 
perpendicular to the scan direction, will cause the new spectra from 
one detector to coincide with those from a different one in the pre- 
vious repeat. 

Each area was observed at least four times, with one repeat offset 
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Table 3. Reference positions used. 
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Figure 2. Representative maps of the RMS noise in K towards IC348 in the 
three CO isotopologues. Top: ^^CO / = 3 ^ 2 noise in 1 kms~^ channels. 
Middle and bottom: ^^CO and C^^O noise in 0.15 km s~^ channels respec- 
tively. 



from the original and the other two basket- weaved variants of the 
original and offset scans. Maps of the resulting root mean square 
(RMS) noise for the three transitions towards IC348 are shown in 
Fig.|2] 



2.2 References and calibrations 

Finding reference positions has proved considerably harder for 
HARP than single-pixel receivers, as all 16 detectors must point 
at line-free regions, requiring a 'blank' piece of sky approximately 
2 arcmin square. Three such positions were used for our Perseus 
survey (see Table|3]). Separate one minute 'stare' observations were 
undertaken at ^^CO / = 3 ^ 2 towards these positions, using ref- 
erences even further from the cloud. All were found to have no 
emission above the noise, re-binned in 1 kms~^ velocity channels, 
in every working detector. 

All the standard telescope observing procedures were fol- 



RA (J2000) 


Dec (J2000) 


Sub-Region used 


(h m s) 


(° ' ") 




03:29:00 


+31:52:30 


NGC1333 


03:48:38 


+31:49:39 


IC348 


03:33:11 


+31:52:02 


L1448andL1455 



lowed. Calibration spectra were taken frequently towards CRL 618 
of the ^^CO or ^^CO / = 3 ^ 2 lines as appropriate. The intensity 
in the reference detector was then compared to previously recorded 
standards and only if they matched within a calibration tolerance, 
were any subsequent observations allowed to continue. All the data 
products and images we present are on the antenna temperature 
[Kutner & Ulich[T98T| 



scale (r^, [Kutner & Ulich|198lf , which can be converted to main 
beam brightness temperature, T^\^, using = T'^/^mb- The effi- 
ciency we use, T7nib = 0.66, was measured during the commission- 
ing of HARP 



3 DATA REDUCTION 

We used the Starlink software collectiorj^for the analysis and re- 
duction of HARP data. First, bad or extremely noisy data were 
flagged in the supplied HARP time- series format. The resultant 
spectra were placed on to a spatial grid using the SMURF reduc- 
tion package ( Jen ness et al.^2008) . Our final data products are sam- 
pled on a 3 arcsec grid, using a 9 arcsec full- width half-maximum 
(FWHM) Gaussian gridding kernel, resulting in an equivalent 
FWHM beam size of 17.7 and 16.8 arcsec for ^^CO/C^^O and 
^^CO respectively. After gridding, a linear baseline was removed 
from each spectrum by fitting line-free portions using the Starlink 
KAPPA applications. The spectra were also re-binned spectrally to 
resolutions of 0.15 km s"^ for ^^CO/C^^O and Ikms"^ for the 
low-resolution ^^CO data. 



3.1 The HARP flatfield 

After the basic data reduction, distinctive stripes were apparent in 
the integrated intensity images. Near the peak of the line intensity, 
rows of high and low pixel values lay parallel to the scan direction, 
implying certain detectors were systematically higher or lower by 
as much as a factor of two. The worst affected data were from the 
^^CO and high-resolution ^^CO scans (see Fig.ji]). A similar pattern 
was clear in every individual scan. 

The origin of such systematic calibration differences between 
detectors is unknown but is likely to be in the intermediate fre- 
quency system. However, for the purposes of this paper we require 
only a pragmatic method to eliminate the striping artefacts from 
molecular cloud data. Therefore, we derive a set of temperature 
conversion factors (TCFs), by which we can multiply the spectra 
from each detector to get them on to a common intensity scale. This 
procedure is similar to the 'flatfield' procedure which accounted 
for the flux conversion factors with SCUBA ^Holland et allT999l ), 
therefore we refer to the analogous flatfield for HARP. 

Ideally, to derive these factors from observations, one would 



Now maintained and developed by the Joint Astronomy Centre, see 



[http: //starlink. jach.hawaii .edu| 
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Figure 3. Integrated intensity images for the different / = 3 ^ 2 transitions 
towards IC348 displaying the striping artefacts from our basic data reduc- 
tion. The spectra are distributed on a 7.3 arcsec grid using nearest-neighbour 
allocation to emphasize detector based variations. From top to bottom are 
C^^O, ^^CO, ^2 CO low-resolution and ^^CO high-resolution maps. The 
^^CO data from a few telescope scan rows appear to have slipped, this is 
a known problem where the telescope does not return to the top of the map 
before commencing integrating. 



observe a small source, providing a fully- sampled map for each de- 
tector, so that the total flux can be compared between them. How- 
ever, we must derive the factors from the scientific observations 
themselves, comparing the intensity detector-by-detector for an en- 
tire map. First, datacubes were produced for each detector individ- 
ually and the spectra integrated over the line and summed over the 
map. We divided the sum from the reference detector (Rl 1, one of 
the central four in the array, see [Buckle et al. 2009 for an expla- 
nation of how the JCMT points and tracks with HARP) by those 
from each of the other working detectors in turn (up to 15 in num- 
ber), to calculate the respective TCF. The reference detector was 
taken as the intensity standard as it is not masked out in any of our 
observations. Additionally, all the calibrations are centred on it, so 
we have a robust determination of its performance. Once the fac- 
tors are known, the output from each detector is multiplied by its 
respective TCF. The procedure worked very effectively for most of 
the data as shown in Fig. |4] and all the subsequent images in this 
paper have been multiplied by TCFs, except the C^^O data towards 
L1455 which had an insufficient signal-to-noise ratio (SNR)and all 
the low-resolution ^^CO datasets, which did not display the arte- 
facts. 

It is perhaps surprising that this procedure worked so effec- 
tively and undoubtedly the nature of the targets contributed to its 
success. The method depends on each detector receiving exactly the 
same amount of emission. Offset map repeats and basket- weaving 
helped each detector to sample a larger proportion of the target 
field. However, it is the uniformity of a molecular cloud's emission 
that means, in general, there is little difference in intensity between 
the different detectors. If there were many compact sources or in- 
tense emission towards the edges of these maps, which are only 
observed by a few detectors, then this technique would be ineffec- 
tive. Other datasets taken in the same period show similar striping 
artefacts, notably those of the JCMT GBS, who have implemented 
our algorithm across a number of their fields. 



4 RESULTS 
4.1 NGC1333 

NGC1333 is a young stellar cluster in the west of Perseus, asso- 
ciated with the reflection nebula of the same name and the dark 
cloud LI 450. It has been widely studied (see [ Walawender et al.| 
|T008 for a review) and is the most active region of star formation 
in the Perseus complex. The stellar cluster is very young (< 1 Myr) 
and highly obscured, containing about 150 stars totalling 79 Mq 
( |Lada, Alves & Lada|[T996l |Wilking et al.""2004V The latest in- 
frared (IR) survey with Spitzer ( J0rgensen et al. 20061 [ Rebull gt al.| 
|2007""Gute rmuth et aL|2008| ) identified 137 objects in the cluster: 
98 pre-main sequence stars and 39 protostars. These protostars are 
correlated with the position of dense molecular material and dust. 
In total, there is approximately 450 of gas in the region f Warin] 
|et al.|1996|). Much of th is gas lies in dense filaments, surrounding 
cavities ( [Lefloch et aL|[T998l [Quillen et ari[2005] ). The YSOs in 
NGC1333 also drive a large number of overlapping outflows, pro- 
viding one of the clearest examples of the self-regulation of star 
formation ( Knee & Sandell|2000| . 

We present maps of the ^^CO and C^^O HARP data in Fig. 
|5] alongside the SCUBA 850 ^tm emission over the observed re- 
gion, originally from observations by 'Sandell & Kneel (| 2001| l but 
later analysed as part of the survey of Hatchell et al. ( 2005|. Qual- 
itatively, the C^^O and SCUBA maps closely resemble each other. 
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Figure 4. Action of the flatfield procedure on ^^CO integrated intensity images towards NGC1333. The panels are from top to bottom and left to right: data 
from a single scan, unaltered from the basic reduction; the same data as in the previous panel with the flatfield applied; data from all the NGC1333 scans, 
unaltered from the basic reduction; and the same data as in the previous panel but with the flatfield applied. 



The central cavity, just north of SVS13 is clear in all the tracers, 
as is the horseshoe of emission enclosing its northern boundary. 
This would suggest that similar material is being explored in both 
the C^^O and SCUBA data. However, there are clear differences 
in the structures as well. IRAS4 and IRAS2 ( [Jennings e t al. 19871, 
are two of the brightest sources in the SCUBA field, yet they have 
much weaker C^^O emission and can scarcely be discerned in the 
^^CO data. The stringof cores to the north of the central cavity, 
(sources 54, 56 and 6qj) have quite weak SCUBA emission but are 
the brightest objects in the C^^O maps. 

4.2 The IC348 molecular ridge 

Our second field is in the vicinity of another young IR cluster 
(-2Myr old), IC348 (see |Herbst|[2008 ^ for a review), in the east 
of Perseus. It comprises several hundred members totalling about 
160 M0 ( [Luhman et al.|2 003l. The low stellar disc fraction and lack 
of outflow activity may indicate that the cluster is coming to the end 
of its star-forming phase ( [Luhman et al.|1998] ). In fact the region 
we refer to as IC348 hereafter (Fig.[6| is really a bright molecular 
ridge some lOarcmin southwest of the cluster, sometimes associ- 
ated with the "Flying Ghost Nebula" (Bou lard et al.|[1995| ). This 
area, in contrast to IC348, is currently undergoing star formation 

^ The source numberings we refer to in the rest of this paper are the cata- 
logue number from |Hatchell et al.[^2007a| 



with many embedded objects and outflows, possibly triggered by 
the nearby cluster (e.g. [Bally et al.|2008] ). The best-known feature 
in the area is HH211, a highly-collimated bipolar outflow driven 
by a Class protostar discovered by'McCaughrean, Rayner & Zin-[ 
[necker[ ( fl994| ), which has been the target of many interferometric 
studies subsequently (e.g. Gueth & Guilloteau[[1999[ [Chandler"&| 
[Richer|2001| ). 

In our C^^O data, depicted in Fig.|6] the close correspondence 
with the SCUBA emission is striking. Again there is not a simple 
scaling between SCUBA flux and C^^O integrated intensity, for 
example the brightest C^^O core is the starless source 17, while the 
driving source of HH21 1 (source 12), the brightest SCUBA core, is 
less prominent in C^^O. The ^^CO data looks very different, as it is 
probably optically thick (thus tracing only outer layers of the cloud) 
and has very bright outflow lobes. Two bipolar outflow structures 
can be seen in the ^^CO integrated intensity and are more obvious 
in the false colour red-green-blue (RGB) image of the same gas. 
The first highly symmetric outflow (driven by source 12) is HH21 1, 
while the second is the more confused north- south flow from source 
13 - kn own also as IC348-MMS ([Eisloffel et al.|2003|) or IC348 - 
SMM2 ( Ta falla, Kumar & B achiller 2006 , Walawender et al.|2006l ). 
The bipolar outflow discovered by Tafalla et al. ( 2006) from IC348- 
SMM3 (source 15), slightly west of IC348-MMS7 is also faintly 
discernible. Two Class I protostars (sources 14 and 101) are very 
red in the RGB image. This may indicate either an outflow where 
only one lobe is visible due to e.g. an inhomogeneous environment 
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Figure 5. Overview o f the data towards NGC1333. Top left: SCUBA 850 i^m emission with contours at 100, 200, 400, 800, 1600 and 3 200mJy beam-| 
( [Hatchell et al.|2005t with a colour-scale as shown in mJybeam ^ . The peak of submillimetre cores identified and classified by H atchell et al. 112007 a) 
are labelled with their [Hatchell et al. | source number and marked with: black and white circle (starless cores), light blue diamonds (Class protostars) or 
purple triangles (Class I protostars). |Sandell & KneeH2001) o riginally presented this data and located many submillimetre sources. Additionally labelled are 
prominent sources in the field, NGC1333- IRAS2A jJennmgse t al. 1987 Blake et al. 1995 Looney, Mundy & Welch 2000), NGC1333-IRAS4A to C (Jennings] 
et al.^1987 Lay, Carlstrom & Hills 1995" Rodriguez, Anglada&Curiel^l999 Looney et al. 2000i, SVS13 (Strom, Vrba & StFom 1976 Haschick et al. 1980 
Snell & Edwards 1981 Grossman et al. 1987 Hirota et al. 200^, HH7-11 ( Herbig 1974) and NGC1333-IRAS7 ( Lightfoot & Glencross 1986 Jennings et al.| 
1987 Cohen, Jones & Hereld 1991 J0rgensen et al. 2007b). Top right: Contours of SCUBA 850 yam emission (as previously) overlaid on the C^^O / = 3 ^ 2 
integrated intensity as displayed on the bottom left. Bottom left: Integrated C^^O / = 3 ^ 2 intensity, /T^dv, from 5 to 10km s~^ with contours from 1 to 
lOKkms"^ in steps of 1 Kkms~^. Bottom right: Integrated ^^CO / = 3 ^ 2 intensity, / T^dv, from 2 to 17kms~^ with contours from 3 to 33 Kkms~^ in 
steps of 3Kkms~^ 
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or that the CO emission comes from the protostars themselves 
and they are moving with respect to the ambient gas. 



4.3 L1448 

The dark Lynds cloud L1448 (see Fig. [7]), the most westerly of 
our targets, is a region dominated by outflow activity. Half a 
dozen or so YSOs reside in its dual-core molecular structure of 
-100 M© (BachiiT er& Cernicharo||l986b| [Wolf-Chase, Barsony| 
|& O' Linger 2000 ). Its outflows are well studied, particularly the 
highly-collimated, symmetric flow originating from L1448-C, one 
of the youngest known at the time (e.g |Bachiller et al.|1990[[Bany^ 
|Lada & Lane 1993 ). Given that the energy in the outflows exceeds 
the gravitational binding energy of the cloud, LI 448 is likely to be 
dispersed by its outflow activity (Wolf-Chase et al.|200 0). 

The HARP C^^O data seems to reflect straightforwardly the 
SCUBA emission in Fig. [7] with the brightest cores at 850 ^im also 
having the greatest C^^O integrated intensity. Most of the SCUBA 
cores also appear as clumps of C^^O emission. The same structure 
is clear in ^^CO as well, even though there is some hint that the 
emission is lined up along the main southeast-northwest outflow 
axis. The outflows themselves dominate the ^^CO data with per- 
haps three to four flows overlapping. L1448C (source 29) drives 
a highly collimated outflow discovered by |Bachiller et al.| ( |1990| ), 
which was one of the highest-velocity and youngest found at the 
time ( |Bachiller et al.|1990|[7995l ). Its northern, blue-shifted lobe in- 
tersects an outflow from a cluster of Class objects: L1448N:A/B 
and L1448NW. Finally, west of this cluster there are perhaps two 
outflows emanating from the vicinity of L1448-IRS2 (sources 30 
and 31). 

4.4 L1455 

Our final field is towards another Lynds cloud, LI 45 5 (Fig. |8]) in 
the south-west of the complex with a mass of some 40-50 
( |Bachiller & Cernicharo| 1986a] ). Although the smallest and faintest 
of our targets, it has some of the most interesting outflow structure. 
The C^^O emission is weak and only detectable towards the hand- 
ful of protostellar and starless cores in the southeast of the field. 
The ^^CO data also mainly pick out the compact sources with as- 
sociated SCUBA emission. One of the brightest ^^CO clumps (also 
detected in C^^O at 3^27"^33^ 30° 12'45'0 has no compact SCUBA 
emission and appears coincident with collimated blue-shifted emis- 
sion from RN015-FIR (L1455-FIR, source 35), a Class I source as- 
sociated with the red reflection nebula of the same name. The ^^CO 
maps are more intriguing, with a prominent northwest- southeast 
CO outflow ( Goldsmith et al. ' 19841 |Levreault| 1988|) some distance 
from the cluster of protostars. Either this is driven by an unknown 
low-luminosity source or one of the known protostars. The latter 
explanation is promoted by |Davis et al.| ( [2008| ) although the struc- 
ture of the outflow might suggest the former is more likely. 



5 DISCUSSION 

5.1 Average region properties 

In Table |4] we list the masses of the regions derived from the CO 
isotopologues and SCUBA 850 \xm data. For the continuum data, 
we follow fVisser, Richer & Chandler] ( [200T] ) in calculating the re- 
gion masses, assuming a constant dust temperature and opacity of 



12 K and K'gso = 0.012 cm^g~^ respectively. The gas masses as- 
sume each spectral line is optically thin, in LTE and its emitting 
gas has an excitation temperature, Tex = 12K. We use a standard 
method to derive the molecular column densities (e.g. |Minchin,| 
[White & Padman|1993| ), which for a CO isotopologue, X, yields: 

A^(X) = 5 X 10l2 Texexp j ^mbdv cm-^, (1) 

where Ttrans is 33.7, 31.8 and 31.6 K for the / = 3 ^ 2 transition of 
^^CO, ^^CO and C^^O respectively and the integral is in Kkm s~^. 
To calculate total masses from the derived column densities we as- 
sume molecular abundances of 10~^, 10~^ and 10~^ for C^^O, 
^^CO and ^^CO respectively relative to neutral hydrogen. A main- 
beam efficiency of = 0.66 was used for all the isotopologues 
as measured during HARP's commissioning. In each map, only po- 
sitions with peak values greater than three times the average noise, 
(orms)» have been included. 

The pattern of the different masses is intriguing and not 
purely what is expected given the different optical depths of 
the CO isotopologues. NGC1333 clearly shows the naively- 
expected trend; as the CO lines become increasingly optically thick 
(C^^O^^^CO^^^CO), only the cloud surface layers are probed, 
the total emission falls and the mass decreases. The SCUBA 
850 ^im emission, although optically thin almost everywhere, is 
less sensitive than the CO observations to large-scale structure 
(> 2 arcmin), due to the spatial chopping in the observing tech- 
nique. Therefore, it potentially misses much of the cloud mass. 
This simple prescription is not followed in any of the other re- 
gions, except perhaps LI 455. However, here the C^^O emission 
is particularly weak so most of the spectra have peaks < 3orms 
and are therefore not included in the mass estimates of Table |4] 
which probably relate to three or four cores. This is a potential rea- 
son why the C^^O mass is so much smaller than that from SCUBA. 
IC348 and L1448 make an interesting contrast: the C^^O and ^^CO 
masses are not very different, with the ^^CO mass marginally larger 
in both cases, implying the ^^CO gas is probably optically thin and 
more widespread than the C^^O. However, LI 448 has a SCUBA 
mass some 3-4 times larger than these gas estimates in contrast 
to IC348 whose SCUBA mass is 2-3 times smaller. C^^O excita- 
tion effects may explain the large C^^O mass in IC348, if there is 
a significant amount of subthermal excitation, the C^^O emission 
for a given quantity of gas would be boosted relative to LTE caus- 
ing an over-estimate of the C^^O mass. In L1448, the small C^^O 
mass is unlikely to be caused by subthermal excitation and proba- 
bly arises from a different gas-to-dust ratio or if we underestimate 
the dust temperature (which in any case will vary from our constant 
assumption of 12 K). In Section [531 we show the dust temperature 
in LI 448 is always higher than the excitation temperature of the 
gas. If we underestimate the dust temperature then we overestimate 
the corresponding dust mass, which may account for some of the 
discrepancy between the tracers. 

Average spectra in the different regions are plotted in Fig. [9] 
All three isotopologues are centred at similar velocities, approxi- 
mately 7.8, 8.8, 4.3 and 5.0 km s-^ for NGC1333, IC348, L1448 
and LI 45 5 respectively. The strength of the lines are in order of 
their abundance as we would expect and these averages show no 
global evidence for multiple components or self-absorption. How- 
ever, individual ^^CO spectra do show self-absorption and/or mul- 
tiple components. In L1448, the ^^CO line is nearly the same 
strength as the ^^CO, implying the ^^CO line is either saturated 
(and therefore the gas has a low physical temperature) or there are 
significant optical depths in the ^^CO gas due to density or tem- 
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Table 4. Region masses from the various HARP / = 3 ^ 2 integrated in- 
tensities and SCUBA 850 lam emission. All spectral-line estimates assume 
LTE, T ^ 1 and Tex = 12 K, while the SCUBA estimate assumes a constant 
dust temperature of 12 K. Only points with peaks > 3{(7rms) have been 
included. 



Usi ng the radiative transfe r equation for an isothermal slab 



Region 




Masses (M©) 






850 lam 


ci^o 


13CO 


12CO 


NGC1333 


250 


412 


301 


13 


IC348 


71 


160 


212 


35 


L1448 


153 


38 


45 


9 


L1455 


45 


7 


35 


12 



perature effects. Finally, even an average ^^CO spectrum contains 
significant high- velocity linewings and there are weaker linewings 
clear in the ^^CO spectra towards L1448 and IC348. 

C^^O is an ideal tracer of the bulk motions of the molecular 
gas, given that it is nearly optically thin throughout and detectable 
over larger areas than denser-gas tracers. We plot C^^O line centres 
across the regions in Fig. [To| The centres are from non-linear least 
squares fitting of Gaussian functions to every spectrum with a peak 
> 3orms (the noise calculated on each individual spectrum), per- 
formed using software from S. Graves. There are many plausible 
explanations for variations in velocity across a cloud region: 

(i) Rotation. A systematic shift in the line velocity along a par- 
ticular direction is often suggested to be a signature of bulk rotation. 
For instance the Serpens cloud core is thought to be undergoing a 
global east to west rotation (e.g. |01mi & Testi|2002) which mani- 
fests itself in a very ordered change in the C^^O J = 3 2 velocity 
observed with HARP by the GBS (Graves et al., in preparation). 
In the Perseus regions there are few such global gradients. |Ho &| 
|Barrett||l980 ) inferred a south to north rotation of NGC1333 from 
NH3 observations, although fWalsh et al.H2007|) (and this wo rk) find 

were bi- 



Ho & Barrett 



no such global gradient in N2H+, suggesting 1 
ased by the red- shifted gas around IRAS 7 (in the centre-left of our 
map). 

(ii) Outflows. On large scales in our fields, CO outflows do not 
seem to correlate with the C^^O centres (although there is some link 
on the scale of individual cores; Curtis & Richer, in preparation). 

(iii) Constituent clumps moving with respect to the bulk cloud. 
Example SCUBA cores seem to enclose regions of very different 
velocity, which implies some of the variations trace clumps of gas 
that move at distinct velocities with respect to the ambient cloud. 

(iv) Gas Flows. In a similar fashion to (iii) the variations in ve- 
locity could reflect larger flows of gas. In gravoturbulent models we 
expect cores to form at stagnation points, where convergent streams 
of material collide (e.g. |Padoan et al.,2001| ). In these data, there are 
large areas of coherent motion but the cores do not seem to occur 
solely where they intersect. 



5.2 ^^CO/C^«0 ratio 

The optical depth of a spectral line provides important evidence as 
to where in the cloud it has been emitted. Optically-thin molecular 
lines can provide a trustworthy measure of the column density and 
conditions throughout the cloud, whereas optically-thick lines satu- 
rate, only tracing the outer layers. From the peak ^^CO/C^^O ratio, 
we can get a measure of the optical depth of the C^^O 7 = 3 ^2 
line in the densest region along the line of sight and thereby test its 
fidelity as a total mass tracer. 



(see e.g. 



Rohlfs & Wilson 



2004| ), we can relate the ^^CO/C^^O in- 
tensity ratio at the velocity of the peak of the C^^O line, R, to the 
optical depths of both species (T13 and Tig respectively), assum- 
ing both transitions are at the same frequency, emanating from the 
same volume of material and have the same T7jnb and beam-filling 
factors: 



R- 



TlC'CO) _ l-exp(-Ti3) 



r;(Ci80) l-exp(-Ti8) 



(2) 



The opacities are expected to be linked through the abundan ce of 
^^CO relative to C^^O, which we denote ^, as T13 = ^Tig (e.g.iMy- 



ers, Linke & Benson[l983| . Where the lines (particularly of ^^CO) 
are optically thin, T ^ 0, and R^ ^ asymptotically. The absolute 
values of the ^^CO and C^^O abundances are not known to high 
precision and may differ across regions or between clouds, due to 
statistical variations across the Galaxy. For instance in photon dom- 
inated regions, R can be very high where there is little self- shielding 
of the cloud from incident radiation, causing the rarer isotopologue 
to be almost completely destroyed (Storz er et al.||2000| . |Wilson| 
|& Rood| ([T997) collated various studies together and tracked the 
changes in relative abundances with distance from the Galactic cen- 
tre. In the local interstellar medium they found [^^C]/[^^C] = 77±7 
and [160]/[1^0] = 560±25, indicating C = 7.3 ±0.7 if there is no 
fractionation. 

In calculations of R, we only include points that have a good 
detection of both the C^^O and ^^CO lines (peak line brightness > 5 
times the RMS spectral noise). We plot R in Fig.[TT] Where the gas 



becomes less dense R tends to C,. In Fig.lllJ outside of the SCUBA 
contours, so presumably at lower space densities, the ratio is ~5-6. 
At the very edges of the measured values in NGC1333 and IC348, R 
reaches ~7-8, the value from | Wilson & Rood| ( |1994| ) and becomes 
larger still. These larger ratios could either be a result of actual de- 
viations in ^ or from noise effects. Most of the largest R are in 
the northeast of NGC1333 and northeast of the main horseshoe of 
SCUBA emission in IC348. These two areas are not representative 
of the two regions. In NGC1333, this is where a bubble of bright 
8.0|j,m Spitzer flux resides (see fig. 1 of I Gutermuth et al.']|2008| ), 
commonly thought to be emission from polycyclic aromatic hydro- 
carbon features excited by UV radiation. The UV source is likely 
to heat up the molecular gas in the region and perhaps change the 
isotopologues' abundances. In the IC348 area, there is a mass of 
blue-shifted ^^CO gas, expanding into a dust cavity ( |Tafalla et af 
2006 and Curtis et al., in preparation). Thus, disregarding these ar- 
eas, R 7 at the edges of the map, which we take to be ^ in the 
following analysis. 

For C, =7, Ti3 = 7ti8 and R varies with t as in Fig. \2 R 
clearly falls towards the centre of the dust filaments and cores, 
where we expect the gas to be more dense. The average R in each 
region is unlikely to relate convincingly the properties of the whole 
area as it is measured over such a small portion of the maps, partic- 
ularly in LI 448 and LI 455. The typical ratios can be measured in 
the majority of the star-forming SCUBA cores, which may prove 
more informative. We measured R at the peak of each core identi- 
fied by |Hatchell et al.| ( [2007a ) and present their averages with those 
in the entire regions for comparison in Tables |5]and[6] 

The typical ratios across the four regions, R^2 — A, trans- 
late into Ti3 p:^ 1.52 - 0.15 and Tig ^ 0.22 - 0.02. This suggests 
that the C^^O line is typically optically thin but the ^^CO can be 
classed neither as optically thin nor thick. The pattern of R suggests 
that on average LI 448 has the highest column density, followed by 
NGC1333, whilst IC348 has the smallest. The value for L1455 is 
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Table 5. ^^CO/C^^O ratios, R, at the peak of the SCUBA cores in the 
[Hatchell et al.|{20"07at catalogue. Gr is the standard deviation of R across 
the cores. 



Region or 
Core Class 


Number 




Or 


NGC1333 


28 


2.8 


0.8 


IC348 


16 


3.0 


0.8 


L1448 


7 


2.2 


0.3 


L1455 


1 


3.6 




Starless 


21 


2.8 


0.6 


Class 


21 


2.6 


0.8 


Class I 


10 


3.1 


1.0 



not representative of the region as it is only measured in the dense 
cores. When just the cores are examined the pattern is the same 
IC348 has the largest R followed by NGC1333 and then L1448. 
As we expect, the cores have larger optical depths than their parent 
clouds on average. However, the reduction in R is not enough to 
render the core C^^O gas optically thick. The smallest core ratio 



Table 6. ^^CO/C^^O ratios, R, for every map pixel in Fig. [IT] Gr is the 
standard deviation of R across the pixels. 



Region 




Or 


Median R 


NGC1333 


3.7 


1.3 


3.5 


IC348 


4.1 


1.8 


3.6 


L1448 


2.8 


0.7 


2.7 


L1455 


3.6 


0.9 


3.8 



is 1.62 in NGC1333, which implies Tig ^1.0 whereas must cores 
have much smaller Tig. The mean core ratios are smaller (and thus 
the column densities are higher) for Class than I protostars, al- 
though there is significant spread in both populations. As collapse 
progresses in protostars, they accrete an increasing fraction of their 
envelopes on to the central object so reducing the column density 
seen towards Class I over cores. The column density seen towards 
starless cores presumably depends on the exact age of the core so 
could be comparable to the protostars' or entirely different. 
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5.3 Temperature 

In this section, we derive Tex for the ^^CO / = 3 ^ 2 hne assum- 
ing it is optically thick. We include ^^CO HARP data for NGC1333 
over the region observed in ^^CO/C^^O (J. Swift, personal commu- 
nication) which are distributed on a coarser grid of 6 arc sec pixels 
using a nearest-neighbour scheme. 

In LTE, Tex is the physical temperature of the gas, although 
this will of course not be the case in general. We have noted typi- 
cal ^^CO optical depths of Tn ^ 1 — 2, which p rovided the abun- 
dance [^2c]/[i3c] = 11 jwilson & Rood|l994| ) should ensure the 



opacity of the ^^CO gas is T12 ^77 — 154. Ratio maps with other 
isotopologues for ^^CO are unlikely to produce accurate T12 val- 
ues as ^^CO traces outer cloud layers and not the volumes where 



are emitted. However, provided the line is optically 
thick, i.e. T ^ 00, and not self-absorbed, we can derive Tex thus 
(see e.g. |Pineda et al.|2008| ): 



rex(''CO) : 



16.59K 



In { 1 + 16.6K/[rmax(^2CO) + 0.036K] } 



(3) 



where Tq = 16.59 K and rmax(^^CO) is the peak ^^CO main beam 
brightness temperature. We plot Tex in Figs. [13] and [14] Tex is typi- 
cally 5-25 K, generally increasing towards the centre of each region 
where it can reach 30 K. 

There are a number of caveats to this analysis that should be 
noted. First, in some respects the most interesting regions of molec- 
ular clouds are the cores in which stars are formed. The tempera- 
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tures of such high-density regions will not be accurately probed by 
^^CO, whose emission is restricted to outer cloud material due to its 
higher opacity. However, we might naturally expect the temperature 
of starless cores (without an internal heating source) and possibly 
protostars to reflect that of their natal environment so the ^^CO tem- 
perature should reflect the pattern of core conditions. Second, ^^CO 
may not be optically thick everywhere, although the typical ^^CO 
optical depths in regions with bright SCUBA emission do suggest 
it will be for most of the areas. Furthermore, we may underestimate 
the true ^^CO peak temperature if the line profile is complicated by 
effects such as infall and self-absorption. 

In L1455, L1448 and IC348, the Tex maps follow the inte- 
grated ^^CO intensities closely. Any strong outflows are promi- 



nent areas of higher Tex, as suggested by other authors ( [Hatchell,] 
[Fuller & Ladd|1999[|Nisrni et al.|2000l ). The dominant feature in 
NGC1333's Tex map is the blue- shifted bubble of gas in the north- 
east corner. We have already mentioned this is an area likely to 
be illuminated by a UV source. By contrast, the most pronounced 
feature in the ^^CO integrated intensity, the bright outflow from 
SVS13, is almost absent. The bubble's Tex is some 40 K higher than 
the maximum in the other regions. In this part of the map, ^^CO is 
not optically thick and the gas probably has a higher physical tem- 
perature - as we noted the ^^CO/C^^O ratio is larger in this region. 
Additionally, it only has single-peaked lines, devoid of some of the 
more complicated profiles seen in NGC1333, suggesting its peak 
temperature is over-estimated compared to the rest of the cloud. 
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Table 7. Tex from ^^CO at the peak of the SCUBA cores in the |Hatchell| 
|et al.|{200 7a ) catalogue by region or source age. 



Region or Core Class 


Number 


(7;x)/K 


(Tr/K 


NGC1333 


29 


25.2 


10.1 


IC348 


17 


17.6 


3.2 


L1448 


7 


10.6 


2.0 


L1455 


5 


12.6 


2.2 


Starless 


23 


18.8 


9.7 


Class 


21 


18.3 


6.2 


Class I 


14 


24.9 


11.2 



This all points to a warm, low-density region with material per- 
haps dispersed by nearby luminous stars, which Hatch ell, Fuller &| 
jRicher ( 2007b ) suggest are the cause of a lack of outflow detections 
in this area. 

It is interesting that the largest temperatures in 
NGC1333 are to the northeast, in the direction of 40 Persei 
(Q3h42m22.6s,+33°5r54.rO, a B0.5 star, part of the Perseus OB 
association, suggested to be triggering the star formation in these 
clouds ( IWalawender et al.|[2004| |Kirk et al.|[2006l ). Indeed, the 
highest temperatures are to the north in IC348, also in the direction 
of 40 Per, although no such gradients are clear in LI 448 and 
L1455, possibly as they are further away. 

On average, LI 448 is the coldest, followed in order of increas- 
ing temperature by L1455 then IC348, with NGC1333 the hottest. 
The histograms of Fig. [14] show two distinct temperature environ- 
ments: one colder in L1455 and L1448 and the other hotter in IC348 
and NGC1333. The latter two regions show a large range of tem- 
peratures. NGC1333 and IC348 are probably heated more by their 
namesake IR clusters and the triggering radiation from 40 Per than 
L1448 and L1455. Protostars may also heat their surroundings to a 
limited extent via radiation and outflows. NGC1333 with its numer- 
ous energetic flows and protostars is most likely to be warmed up 
in this way. LI 448 has a high protostellar fraction as well and pow- 
erful outflows but less ambient gas with which to interact - most of 
the outflows have broken out of the region. 

As we have already mentioned, temperatures derived from 
^^CO are unlikely to be best probes of conditions inside dense 
cores. Nevertheless we might expect the core temperatures to be 
affected to some extent by the temperature in the bulk of the gas. 
From Table |7] it is clear the temperatures towards the cores are 
larger on average than within their parent region as a whole but 
exhibit the same trends. When comparing cores of different ages 
it is worth noting that we expect fundamentally different relations 
between Tex and the true temperatures for starless and protostellar 
cores, Tcore- Towards the centres of starless cores, temperatures fall 
whereas for protostars temperatures rise. Thus, for starless cores 
Tex > Tcore and for protostars Tex < Tiore- In Table [7] the environ- 
ments of Class I protostars appear slightly warmer than both star- 
less cores and Class protostars on average. As a source ages we 
expect its core temperature to increase, resulting from the increase 
in brightness of its central object. These data would exhibit this 
trend except the Class and starless cores have similar tempera- 
tures, probably illustrating that the ^^CO line probes an environ- 
ment outside the high-density core, where Tex > Tiore for the star- 
less cores and Tex < Tiore for the protostars. 

We can make an interesting comparison between this gas tem- 
perature and that of the dust. At high densities, n>2x 10^ cm~^, 
conditions likely to be probed in CO with HARP, the gas and dust 



temperatures are couple d ([Galli, Walmsley & Gongalves||2002| ). 
Recently, [Schnee et al.| ( |2008| ) computed the dust temperature at 
40 arcsec resolution across Perseus using Spitzer and IRAS mea- 
surements. We plot Tex degraded to 40 arcsec resolution versus 
[Schnee et al.[ s dust temperature in Fig. [15] The dust temperature, 
7b ^ was calculated from Spitzer 70 and 160 ^im fluxes, which is 
dominated by warm dust along the line of sight rather than the 
denser, colder material in say starless cores and is more likely to 
match the material traced in ^^CO. There is actually little correla- 
tion between the two temperatures. Typically Td is slightly lower 
than Tex: 12-20 K compared to 8-30 K. In NGC1333 and IC348, 
Tex is mostly larger than although there is a significant propor- 
tion of pixels where the opposite is true. However, in LI 448 and 
LI 455, nearly everywhere Td > Tex- This may be explained if the 
^^CO gas is not tracing the same areas as the dust emission, with 
the ^^CO emitted in a colder layer whilst the dust comes from a hot- 
ter inner cloud region. Additionally, if the gas is not optically thick, 
then we would over-estimate the ^^CO excitation temperature and 
Tex > Tb- 

The ^^CO excitation temperature may be correlated with the 
cloud column density. As the rotational transition is collisionally 
excited, on moving to denser and higher extinction portions of the 
cloud with densities above the critical density, the gas Tex should in- 
crease. We use the integrated C^^O/ = 3^2 intensity as a proxy 
for the visual extinction since typically there is a linear relation be- 
tween them. Using C^^O/=1^0 data across the Perseus molec- 
ular cloud Pineda et al.| ( |2008| ) measured: 

Av/mag = (2.4±0.1)f ^-^^^^^^^ +{2.9 ±0.9). 



V Kkms 



(4) 
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we plot the C^^O integrated intensity, / dv versus 



In Fig. 

Tex, with lines denoting the expected variation for various column 
densities of C^^O assuming LTE. No single line is a good fit for 
all the points as in each region there is a range of column densities. 
However, a moderate range of column densities will span the whole 
parameter space. At low temperatures with a good deal of scatter 
the integrated intensity increases linearly with Tex- LI 448 again re- 
quires the highest column densities with IC348 and NGC1333 hav- 
ing a large spread in conditions. 



6 SUMMARY 

This paper presents the technical details and preliminary analy- 
sis of a large-scale survey of the kinematics of molecular gas in 
the Perseus molecular cloud. Observations of the 7 = 3 ^ 2 rota- 
tional transitions of ^^CO, ^^CO and C^^O in over 600arcmin^ of 
NGC1333, IC348, L1448 and L1455 were undertaken with HARP 
on the JCMT. We introduce a new 'flatfield' procedure to account 
for striping artefacts in HARP scan maps of molecular clouds, ap- 
parently resulting from differential performance across the detec- 
tors and/or their samplers. The data from each working detector 
is multiplied by a conversion factor which scales its intensity to 
match the nominated reference detector. The factors are computed 
from the scientific observations themselves by calculating the total 
intensity received by each detector across the whole map. 

We compare integrated intensity maps of the three tracers to 
SCUBA 850 ^im emission (Ha tchell et al.|2005| and the position of 
protostars and starless cores IHatchelTet al.||200"7a| in each field. 
There is a striking similarity between the SCUBA maps and C^^O 
emission, hinting that similar densities of material are traced with 
the gas and dust strongly coupled. However, the detailed structure 
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of individual star-forming cores is not always so simple, for ex- 
ample very bright SCUBA cores sometimes have only weak C^^O 
emission and vice versa. Many outflows are obvious in the ^^CO 
data, which we will examine in detail in a subsequent study (Curtis 
et al. in preparation). The ^^CO maps are somewhat intermediate 
between the other two isotopologues and average spectra across the 
maps do show weak linewings from outflows. Intriguingly, masses 
derived from the HARP and SCUBA data exhibit different trends 
across the four regions, emphasizing the need to examine the de- 
tailed excitation conditions across regions rather than simple con- 
stant assumptions. 

From the ^^CO/C^^O integrated line ratio, R, we explore vari- 
ations in the gas opacity and estimate the two species' relative 
abundance ([^^CO]/[C^^O] - 7). Across most of the regions the 
C^^O gas is optically thin (Tig = 0.02 — 0.22) and therefore is a re- 
liable total mass tracer. Indeed, inside the denser parts of the clouds, 
the star-forming cores, the opacity does not increase beyond a max- 
imum, Ti8 = 0.9. When we also consider that the critical density of 
the transition is some lO^cm"^, we expect to probe an interme- 
diate region between a dense core and its envelope with the C^^O 
line. The ^^CO optical depths are neither in the optically thin nor 
thick regimes (T13 is typically 0.15-1.52). Class protostars have 
smaller ratios than Class Is as expected as more material accretes 
on to the central object over time. 

If we assume the ^^CO line is optically thick, an estimate of 
the excitation temperature can be gathered from the peak line tem- 
perature, assuming LTE (in which case Tex is the physical gas tem- 
perature). In general we derive temperature of 5-25 K, increasing 
towards the centres of the individual regions and in outflow lobes. 
An area in the northwest of NGC1333 has temperatures of over 
45 K, probably as it is heated by a UV source and/or the region has 
lower density so is not optically thick. IC348 and NGC1333 are on 
average much warmer than LI 448 and LI 455. This is partly as the 
averages for NGC1333 are skewed upwards because of the warm 
gas in the north and perhaps since NGC1333 and IC348 are closer 
to 40 Per (thought to be triggering star formation in Perseus). There 
is little correlation between Tgx and the dust temperature, To, de- 
rived from Spitzer observations ( |Schnee et al.|2008] ). Typically To 
is slightly lower than Tq^\ 12-20 K compared to 8-30 K. However, 
in LI 448 and LI 445, typically Tex < 7b, this may imply we are 
tracing completely different regions with the ^^CO gas and warm 
dust or that the gas is not optically thick, so we are over-estimating 

Tex- 

This work demonstrates the utility of HARP for large-scale 
surveys of gas kinematics in nearby molecular clouds. The J = 3^ 
2 transitions of CO and its isotopologues are powerful probes of the 
conditions of star formation when used in combination, examining 
moderately high densities comparable to the dust densities seen by 
SCUBA. 
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Figure 6. Overview of the data towards IC348 as for NGC1333 in Fig.|5] 
Top: SCUBA 850 lam emission. Prominent sources in the field are labelled: 
HH211 jMcCaughrean et al. 1994 Gueth & Guilloteau 1999 Chandler 
Richer 2001 1, IC348-MMS (Eisloffel et al. 2003 Tafalla et al. 2006, 
^ Walawender et al. 2006 1 and IC348-SMM3 (Tafalla et al. 2006 1. Mid- 
dle: Integrated C^^O / = 3 ^ 2 intensity from 7 to lOkms"^ with con- 
tours from 0.5 to 3.0Kkms"^ in 0.5Kkms~^ steps. Bottom: Integrated 
^^CO / = 3 ^ 2 intensity from 5 to 12kms~^ with contours from 2 to 
18Kkms~^ in2Kkms"^ steps. 
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Figure 6 - continued Top: Contours of SCUBA 850 lam emission (levels 
as previously) overlaid on the C^^O / = 3 ^ 2 integrated intensity as be- 
fore. Middle: Integrated ^^CO intensity from to 20km s~^ with contours 
from 15 to 75Kkms~^ in steps of 5Kkms~^. Bottom: Red-green-blue 
colour composite image of the mean ^^CO data value in various velocity 
ranges: from 3.5 to 5.5 km s^^ (blue), 7.0 to 9.0kms~^ (green) and 12.0 to 
14.0kms-i (red). 
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Figure 7. Overview of the data towards LI 448 as for IC348 in Fig. |6] 
Top: SCUBA 850 lam emission. Prominent sources in the field are labelled: 
L1448NW, L1448:A/B (C uriel et al.|1990|[l999^ , L1448C ^Bachiller etal.| 
[T990IIT995) and L1448-IRS2 jWolf-Chase et al.|2000) . 
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Figure 7 - continued Top: Contours of SCUBA 850 i^m emission (levels as 
previously) overlaid on the C^^O / = 3 ^ 2 integrated intensity as before. 
Middle: Integrated ^^CO intensity from —25 to 25kms~^ with contours 
from 10 to 90Kkms~^ in steps of lOKkms"^. Bottom: Ranges: from —25 
to Okms~^ (blue), 2 to 6kms~^ (green) and 7 to 14kms~^ (red). 
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Figure 8. Overview of the data towards L1455 as for IC348 in Fig. [6] Top: 
SCUBA 850 \im emission. Additionally labelled is L1455-FIR ( [Davis et al.| 
[1997). 
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Figure 8 - continued Top: Contours of SCUBA 850 ]xm emission (levels as 
previously) overlaid on the C^^O / = 3 ^ 2 integrated intensity as before. 
Middle: Integrated ^^CO intensity from —5 to 15 km s~^ with contours from 
15 to 65Kkms~^ in steps of 5Kkms~^. Bottom: Ranges: from —5.0 to 
3.5kms-i (blue), 5.0 to 6.5 km s'^ (green) and 8.0 to 14.0 km s'^ (red). 
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Figure 9. Average cloud spectra for the different regions. Each plot con- 
tains the averages of every spectrum in the HARP / = 3 ^ 2 cubes: ^^CO 
(dashed), ^^CO (dot-dashed) and C^^O (solid). 
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Figure 12. Variation of R with t for ^ = 7. 
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Figure 15. Points mark pixels in the various regions' maps: NGC1333 
(blue), IC348 (red), L1448 (green) and L1455 (orange). Top: Tex derived 
from the peak ^^CO versus the dust temperature derived from Spitzer 
MIPS data at 70 and 160 ^im with 40 arcsec resolution ( Schnee et al.pOOS) . 
The line marks where both temperatures are equal. Bottom: C^^O integrated 
intensity versus ^^CO derived T^x. The lines denoted the anticipated depen- 
dence for constant column densities of C^^O as labelled. 



